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Abstract 

We construct the gravitational energy-momentum of the Bondi- 
Sachs space-time, in the famework of the teleparallel equivalent of 
general relativity (TEGR). The Bondi-Sachs line element describes 
gravitational radiation in the asymptotic region of the space-time, 
and is determined by the mass aspect and by two functions, c and d, 
that yield the news functions, which are interpreted as the radiating 
degrees of freedom of the gravitational field. The standard expression 
for the Bondi-Sachs energy-momentum is constructed in terms of the 
mass aspect only. The expression that we obtain in the context of 
the TEGR is given by the standard expression, which represents the 
gravitational energy of the source, plus a new term that is determined 
by the two functions c and d. We interpret this new term as the energy 
of gravitational radiation. 
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1 Introduction 


One of the most interesting consequences of general relativity is the descrip¬ 
tion of gravitational radiation generated by isolated astrophysical configura¬ 
tions, that loose energy in the form of radiation. The space-time around these 
configurations is not strictly asymptotically flat because although the metric 
tensor components fall off in the expected way as g, w = r ] fiu + h^l/r) in the 
asymptotic limit r —> oo, the time derivative of h^ v is of order 1 jr at spacclikc 
infinity in Cartesian coordinates. The first significant work on this issue is 
due to Bondi and collaborators [I], who established the coordinates and no¬ 
tation that are currently employed in the analysis of gravitational radiation. 
The line element obtained by Bondi and collaborators is constructed out of 
the functions M(u, 9) and c(u, 9), where u is the retarded time (u = t — r), 
and r and 9 are spherical coordinates. These functions do not depend on the 
spherical coordinate 0 because the line element is axially symmetric. The 
function M is called the mass aspect, and the time derivative doc = dc/du 
is identified as the first news function. 

Bondi’s line element was subsequently generalised by Sachs [2], who aban¬ 
doned the axial symmetry and obtained the most general metric tensor that 
describes gravitational radiation at spacclikc and null infinities. In the work 
by Sachs there also appears the function d, that yields the second news func¬ 
tion dod. The two news functions are interpreted as the radiating degrees 
of freedom of the gravitational field. The metric tensor obtained by Sachs 
yields what is presently known as the Bondi-Sachs space-time. The mathe¬ 
matical expression of the Bondi-Sachs metric tensor is rather intricate. There 
are very good review articles that clarify the several aspects of the subject, 
and also explain the emergence of the Bondi-Sachs energy-momentum vector 

Mia El HUE]. 

The Bondi-Sachs energy-momentum is constructed out of the mass as¬ 
pect M(u,9,<p ) only, i.e., it does not depend on the functions c(w, 0, 0) and 
d(u,9,(j >) (see Eq. (4.4) of Ref. [6J). To some extent, it is intriguing that 
the total energy and momentum do not depend on the these functions. Re¬ 
cently, the total Bondi-Sachs energy-momentum has been compared to the 
ADM expression [9] of the total gravitational energy-momentum at spacclikc 
infinity of the Bondi-Sachs space-time. The metric tensor of the latter has 
been rewritten in the ordinary (t, r , 9, 0) coordinates by parametrizing the 
spacelike hypersurfaces by the standard time t. It has been found mmm 
that the resulting expression for the total ADM energy-momentum depends 
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on the mass aspect, as expected, but also depends on the functions c and 
d. Although the total ADM energy-momentum is strictly constructed for 
asymptotically flat space-times, the deviation of the resulting expression from 
the standard Bondi-Sachs energy-momentum - and its dependence on the 
functions c and d - is a very interesting result. 

In this article we obtain the expression of the total gravitational energy- 
momentum of the Bondi-Sachs space-time in the context of the teleparallel 
equivalent of general relativity (TEGR) [13]. The latter is an alternative 
geometrical formulation of general relativity based on the tetrad held. The 
TEGR provides a natural geometrical setting for consistent definitions of 
energy, momentum and angular momentum of the gravitational held. The 
definitions are given by surface integrals, they satisfy conservation equations, 
the algebra of the Poincare group, and arise from well defined densities [13]. 
These definitions are possible because of the structure of the held equations 
and the covariance of the tetrad helds under global SO(3,l) transformations. 
The expressions for the gravitational energy-momentum and angular momen¬ 
tum are covariant under global SO(3,1) transformations. In special relativity 
the energy-momentum and angular momentum of localised material systems 
are frame dependent, and so they are in the TEGR, since the presence of 
gravitational helds (in the Newtonian approximation, for instance) do not 
modify this situation. 

We find that the total energy-momentum of the Bondi-Sachs space-time 
is given by the standard integral of the mass aspect, plus a new term that 
contains the time derivative of the functions c and d. This new term is 
interpreted as the energy of gravitational radiation. Attempts have been 
made in the past to arrive at such quantity, but failed because there was no 
guarantee that the integrals would be convergent [3j. The expression that we 
obtain in the present geometrical framework is finite. We test our expression 
by using a simple expression for the news function suggested long time go in 
the literature. 

In section II we briefly present the geometrical framework of the TEGR, 
the definition of the energy-momentum 4-vector, and the conservation equa¬ 
tions. The Bondi-Sachs line element is presented in section III. In this section 
we display the asymptotic expansion of the metric tensor components and 
also construct the set of tetrad helds that will be used in the calculations. 
Section IV contains all relevant expressions and steps that yield the gravita¬ 
tional energy-momentum of the Bondi-Sachs space-time. The application to 
the news function suggested by Papapetrou [2D], Halliday and Janis [21 j and 
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Hobill [22] is given in section V. Finally in section VI we present the final 
remarks. 

Notation: space-time indices n,u,... and S0(3,l) (Lorentz) indices a,b,... 
run from 0 to 3. Time and space indices are indicated according to // = 
0, i, a = (0), («). The tetrad fields are represented by e%, and the torsion 
tensor by T atll/ = d^e a u — d u e af i- The flat space-time metric tensor raises and 
lowers tetrad indices, and is fixed by rj ab = e a ^e bu g 111 ' = (—1,+1,+1,+1). The 
frame components are given by the inverse tetrads {e 0 M }. The determinant 
of the tetrad field is written as e = det(e%). 

The torsion tensor T aiiv is sometimes related to the object of anholonomity 
fi A ny via n,, = e a x T a It is important to note that we assume that the 
space-time geometry is defined by the tetrad fields only, and thus the only 
possible non-trivial definition for the torsion tensor is given by T a tiv . This 
tensor is related to the antisymmetric part of the Weitzenbock connection 
T x v = e aX d^e aV) which determines the Weitzenbock space-time and the dis¬ 
tant parallelism of vector fields. 

2 A brief review of the TEGR 

The teleparallel equivalent of general relativity is a theory for the gravita¬ 
tional field based on the tetrad field. The dynamics of the gravitational field 
in the TEGR is exactly the same as in the standard metric formulation of 
general relativity. The tetrad fields have 16 independent components, and 
the extra six components, compared to the ten components of the metric 
tensor, allow the establishment of additional geometric structures. In the 
geometrical framework determined by the tetrad fields, one may dispose of 
the concepts of both the Riemannian and Weitzenbock geometries. The 
equivalence of the TEGR with Einstein’s general relativity is established by 
means of an identity between the scalar curvature R(e), constructed out of 
the tetrad fields, and a combination of quadratic terms of the torsion tensor, 

eR(e) = -e(-T abc T abc + l -T abc T bac - T a T a ) + 2 d^eT *). (1) 

The formulation of Einstein’s general relativity in the context of the telepar¬ 
allel geometry is discussed in several references, see m M ESI ESI EH M- 
The Lagrangian density of the TEGR is given by the combination of the 
quadratic terms on the right hand side of Eq. (1), 
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L 


= -ke(-T abc T abc + l -T ahc T bac - T a T a ) - L M 
= -keYR bc T abc - L m , (2) 


where k = c 3 /167tG, T a = T b ba , T abc = e b ^e c v T atlVl and 


Y^obc _^-pabc _|_ rj-i b ac rpcab ^ _|_ — (fj aC 'T b Tj ab T C ) 


(3) 


Lm represents the Lagrangian density of the matter fields. The field equa¬ 
tions derived from (2) are equivalent to Einstein’s equations. They read 


e a xe b M^ bXu ) ~ e(Z b " a T bvil - 1 e afl T bcd E bcd ) = ±-eT afl , (4) 

where 5L M /8e a ^ = eT ap . It is possible to show that the left hand side of the 
equation above may be rewritten as \e R af i(e ) — \e a ^R(e) , which proves 
the equivalence of the present formulation with the standard metric theory. 
Equation (4) may be rewritten in a simplified form as 


d u (e'E aXu ) = —e e a ^ + T x ») , (5) 

where T XfJl = e a A T a/i , and t x ^ is defined by 

t Xfi = k(AYj bcX T bc ^ - g x ^ bcd T bcd ). (6) 

In view of the antisymmetry property £ a/ “ / = — it follows that 


d. 


ee%(f v + T A/i ) 


= 0 . 


The equation above yields the continuity (or balance) equation, 

j t J^d 3 xe e%(t° M + T° M ) = ~f s dSj [e + T”) 

We identify f A/i as the gravitational energy-momentum tensor pH], and 


(7) 


( 8 ) 


P a = j d 3 x e e\(t 0fi + T° M ), (9) 

as the total energy-momentum contained within a volume V of the three- 
dimensional space. In view of (5), Eq. (9) may be written as 
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( 10 ) 


P a = -J d 3 xd,n ai , 

where n aj = —4/ceE a0j , which is the momentum canonically conjugated to 
e a j [13j [16] HE] • The quantity djU aj is a well defined space-time scalar density, 
that transforms as a vector under the global SO(3,l) group. The expression 
above may be transformed into a surface integral, and is the definition for 
the gravitational energy-momentum discussed in Refs. mmm, obtained 
in the framework of the Hamiltonian vacuum held equations. Note that Eq. 
(8) is a true energy-momentum conservation equation. 

The emergence of a non-trivial total divergence is a feature of theories 
with torsion. The integration of this total divergence yields a surface integral. 
If we consider the a = (0) component of Eq. (10) in spherical coordinates 
and a spacelike surface S determined by r = constant, we have 

P (0) = — = ~<f dSj n (0)i =4 k<f dddcp £ (0)01 . (11) 

C J s J s 

Adopting asymptotic boundary conditions for the tetrad fields, we find [19] 
that in the limit S —> oo the resulting expression is precisely the surface 
integral at spacelike infinity that defines the ADM energy [9]. This fact is 
a indication that Eq. ([9]) does indeed represent the gravitational energy- 
momentum vector. But note that there is no restriction regarding the appli¬ 
cability of definition (fTOj) . The latter may be applied to arbitrary space-times, 
with arbitrary boundary conditions. 


3 The Bondi-Sachs space-time and the tetrad 
fields 

The Bondi-Sachs metric tensor describes gravitational radiation at null and 
spatial infinities, and in both asymptotic limits the the metric tensor () jW 
approaches the flat space-time metric tensor rj liu , but time derivatives of 
fall off as 1/r. As a consequence, gravitational waves may in principle 
be detected at spatial or null infinities. The line element of the Bondi-Sachs 
space-time in spherical ( u , r, 9, 0) coordinates, where u = t — r is the retarded 
time, is constructed out of the functions M(u,9,(p), c(u,6,(p ) and d(u,9,(ft). 
It is given by 
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ds 2 = g 00 du 2 + g 22 dO 2 + g 33 d(p 2 

+2g 0 i du dr + 2g 02 du dO + 2 g 0 3 du dtp + 2g 23 dd dp , (12) 

where 

—e 2/3 — r 2 (e 27 I/ 2 cosh 25 + e _27 W 2 cosh 25 + 2UW sinh 25 ), 
r 

-e 2 ^, 

—r 2 (e 27 U cosh 25 + W sinh 25 ), 

—r 2 sm.0(e~ 2l W cosh 25 + U sinh 25), 
r 2 e 27 cosh 25 , 
r 2 e ~ 27 cosh 25 sin 2 0 , 

r 2 sinh 25 sin 9 . (13) 

We adopt the usual convention (u, r, 9, p) = (x°, x 1 , x 2 , x 3 ). The functions (3, 
7 , 5, U and W in the equations above are given only in asymptotic form, in 
powers of 1 jr. 

In this article we will present the asymptotic expansions of all held quan¬ 
tities that are effectively needed in the calculations, i.e., we will dispense 
with the powers of 1/r of the held quantities that do not contribute to the 
calculations. Thus, the asymptotic form of the functions above are 


fl'oo — 

fi'oi = 
#02 = 
#03 = 

922 = 

933 = 

923 = 


where 


V ~ 

/3 — 

ry 

5 ~ 
U ~ 
W — 


—r + 2 M , 

c 2 + d 2 

4 r 2 ’ 
c 

r ’ 
d 

~T ^ 

K u i e A) 

r 2 ’ 

r 2 ’ 


l = d 2 c + 2c cot 9 + d 3 d esc 9 , 


( 14 ) 
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I = d 2 d + 2 d cot 9 — d 3 c esc 6. 

In the limit r —>■ oo, the asymptotic form of the functions above yield 


9oo 

r^j 

1 

+ 

to 

9oi 

r^j 

r 

c 2 + d 2 
+ 2r 2 ’ 

902 


l -\ — (2 cl + 2 dl — p ), 
r 

903 


Ism 9 H—(—2 cl + 2 dl — p) sin 9 , 
r 

922 


r 2 + 2 cr + 2(c 2 + d 2 ), 

933 


[r 2 — 2 cr + 2(c 2 + d 2 )} sin 2 6 , 

923 

r^j 

4 d^ 

2 dr sin 6 -\ -sin 6 . 

3 r 


The functions p and p are defined in Refs. mm They depend on functions 
that are not defined above, but since they will not contribute to the final 
expressions, we will not present their definitions here. 

The expressions of the contravariant components of the metric tensor are 
calculated by means of the standard procedure out of Eqs. (fT3j) (not out of 
Eqs. (fT5jl ). The inverse components are given by g l ' u = (—1) M+I/ (1 /g)M )iV , 
where g = —^7oi (^ 722^733 — . 92 . 3 ) is the determinant of the metric tensor, and 
v is the co-factor of the pv component (we have taken into account all 
necessary powers of 1/r of the functions given in Eq. (fTTjl h We find 


,00 


.01 


.11 


,12 


.13 


.22 


g 02 — g 03 — o, 

c 2 + d 2 

-1 2r 2 ’ 

2 M 


l 

^2 ’ 
l sin 6 


1 
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,33 


1 


9 



9 


,23 


( 16 ) 


Now we turn to the construction of the tetrad fields. The inverse tetrads 
e a ^ determine the frame adapted to a particular class of observers in space- 
time. Let the curve x^{r) represent the timelike worldline C of an observer 
in space-time, where r is the proper time of the observer. The velocity 
of the observer along C is given by u M = dx^/dr. A frame adapted to 
this observer is constructed by identifying the timelike component of the 
frame e(o) M with the velocity of the observer: e(o) ^ = u^{r). The three 
other components of the frame, ep) are orthogonal to e(o) M , and may be 
oriented in the three-dimensional space according to the symmetry of the 
physical system. A static observer in space-time is defined by the condition 
— (u 0 , 0,0,0). Thus, a frame adapted to a static observer in space-time 
must satisfy the conditions e( 0 )*(t,a; fe ) = (0,0,0). It is easy to verify, by 
means of a coordinate transformation, that in terms of the retarded time u 
we also have e(o) l (u, x k ) = (0, 0, 0). The Bondi-Sachs space-time is not axially 
symmetric, and therefore there are no distinguished directions at spacelikc 
infinity. Since we will evaluate surface integrals at spacelike infinity, i.e., 
we will be interested only in total quantities (we will also integrate over a 
surface S determined by r = constant, for r finite but sufficiently large), 
any set of tetrad fields that satisfy the asymptotic expansion e a/i ~ r] ail + 
(l/2)h ail (l/r) in Cartesian coordinates when r —>■ oo, and that satisfy the 
conditions e^ l (u,r,9,(p) = 0, will serve our purposes. For such a frame, 
e (i) ^ e ( 2 ) M and e( 3 ) ^ will define the usual unit frame vectors in the x, y and 
0 directions, respectively, in the limit r —» oo, provided e a ^ is constructed 
in Cartesian coordinates. If we restrict the Bondi-Sachs metric to the Bondi 
metric tensor (by making d = 0 = /), then the latter is axially symmetric 
and e( 3 ) M will define the unit vector in the z direction at spacelike infinity. 

It is not straightforward to construct a simple set of tetrad fields that 
yields Eq. (1T2]) . and that satisfy the conditions e( 0 ) 1 = 0. Note that e( 0 ) * = 0 
implies o — 0. One such set of tetrad fields that satisfy these requirements, 
and acquires the asymptotic form e afJ , — r/a^ + (1/2 )h a/1 at spacelikc infinity, 
is given in (u, r, 9, 0) coordinates by 


e (o 


(-A -A, —F, -G ), 




e (i)/i = (0, Bi sin 9 cos 0 + B 2 cos 9 cos 0 — B 3 sin 9 sin 0 , 

C 3 r cos 9 cos 0 — C 2 sin 9 sin 0 , 

—Dr sin 9 sin 0) , 

e ( 2 )/i = (0, Bi sin 9 sin 0 + B 2 cos 9 sin 0 + B 3 sin 9 cos 0, 

C\r cos 9 sin 0 — C 2 sin 9 cos 0 , 

Dr sin 9 cos 0), 

e ( 3 )/i = (0, Bi cos 6* — B 2 sin 9, — C 3 r cos 9, 0). (17) 

The quantities A, B 3 , B 2 , B 3 , Ci, C 2 , D , E, F, G are determined by requir¬ 
ing that e afl yields the metric tensor components (fT5|) according to e afJ ,ebu r ] ab = 
g tw . The determination of exact form of these quantities in terms of the met¬ 
ric tensor components (fl3|i is very complicated and useless for our purposes. 
We will need the components of the torsion tensor only in the asymptotic 
limit r —» 00 . These quantities must satisfy the following equations, 


- A 2 

— 9oo 


-AE 

= 9 01 


-AF 

= 902 


-AG 

= 903 


—E 2 + B\ + B 2 + Bl sin 2 9 

= 9u 

= 0 

EF + B 2 {Cir) + B 3 C 2 sin 2 9 

— 9i 2 

= 0 

-EG + B 3 {Dr) sin 2 9 

= .913 

= 0 

—F 2 + (Cir) 2 + C 2 sin 2 9 

= 922 


-G 2 + (Dr) 2 sin 2 9 

— 933 


-FG + C 2 (Dr) sin 2 9 

— 923 



As we mentioned earlier, in the asymptotic expansion of the field quan¬ 
tities we will display the terms only up to the power of 1/r that is actually 
needed in the calculations, taking care that we do not neglect any relevant 
term up to (1/r) 2 . Terms of order (l/r) n , with n > 3, do not contribute to 
the final, total expressions. We find 
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E 


F ~ 
G ~ 
B\ ~ 
B 2 — 
B 3 ~ 

- 

G 2 - 
D — 


M 

n— > 

r 

—l - (2 cl + 2 dl + Ml — p ), 

r 

l H— (—2 cl + 2dl + Ml — p) 
r 


■ sin 6* 
M 


1 + 

r 

— - — \{2Ml + cl — p), 


1 

sin 0 


—I—-(2 Ml — cl + 2c?/ — p) 


i c 1 

1 + - + ^T 


r / 2 


+ c 2 — cr 


-- 2c? +-(// + 2cc?) 

sm y L r 

1 c 1 (P 2 s 

1-1—+ c -\~ d* 


(19) 


The expressions above completely fix the set of tetrad fields given by Eq. 
(ITT]) . It must be noted that the sin 6 in the denominator of B 3 and C 2 in 
the expressions above does not imply a divergence of e a/J on the z axis of the 
coordinate system, when 6 = 0. In Eq. (1171) these quantities are multiplied 
by sin#. It can be shown that in ( u,x,y,z ) or in ( t,x,y,z ) coordinates, the 
set of tetrad fields given by Eq. (ITT)) is everywhere smooth in the three- 
dimensional space. 


4 The gravitational energy-momentum 

4.1 Gravitational energy 

The gravitational energy contained within a two-dimensional spacelike sur¬ 
face S (defined by r = constant) in an arbitrary space-time, in spherical 
coordinates, is obtained by evaluating the quantity E^ 0 -* 01 , as indicated in 
Eq. (HID . The simplification of this quantity is crucial to arrive at the final 
result. Taking into account Eq. (J3J), and the fact that g m = g 02 = g 0i = 0, 
we find 
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£( 0)01 


_ g ( 0 ) ^£^01 

= e<°> Is 001 + e (°> xE 101 + e® 2 E 201 + e (°> 3 E 301 


= A 
+ E 
+ F 
+ G 


'I 

A 

'I 

A 

'I 

A 

'I 

A 



| jiOOl 
| /j-iOll 
| rj-i02 1 
j rj~i 031 



+ ls m T° 

+ \(g n T’>-g’>'T l ) 
+ -g 21 T° 

2 y J 

+ -g 31 T° . 

2 y J 


( 20 ) 


where = T Q Q M . 

In Eq. (fTTTh E^ 01 is multiplied by the determinant e, whose asymptotic 
expression is e ~ r 2 sin 6. In order to obtain a non-vanishing value for the 
gravitational energy when the surface S of integration approaches the limit 
r —)■ oo, we must select the terms of E*^ 01 that are of the order 1/r 2 . Terms 
of the order 1 /r n , with n > 3, do not contribute to the final expression. 

By expanding Eq. (12UD in terms of the covariant torsion tensor T ajllJ , we 
observe that all terms will be of the type g^ a g° 13 p l7 T a ^ 7 . In view of Eq. (fTEjh 
we see that all products of the type gf ia g 0l3 g 1 'y that arise in the expansion 
of Eq. (l20l) are at least of the order 1/r 2 . Products of the order l/r n , with 
n > 3, yield a vanishing contribution in the limit r —> cxo. Therefore we keep 
only the products that are of the order 1/r 2 , and consider the contributions 
from T afll/ whose values at spacelikc infinity are of the order l/r°. It is 
important to mention that there does not arise any divergent term (of order 
0(r) or higher) in the expansion of eE 1 - 0 ^ 01 . The hnal expression turns out to 
be finite. After long calculations, we conclude that the non-vanishing value 
of E(°) 01 are simplified to 


S(°) 01 = a E°°i + E Eioi + F E 2 °i + G E 3 °i 
= (g 01 g 01 g 22 T 212 + g m g m g 33 T u ^) 

+^E (5f 01 5f 01 5f 22 T 2 02 + fi' 01 fi' 01 fi' 33 r 3 03) 
+\Fg 01 g 01 g 22 (T 012 -T 201 -T W2 ) 
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(21) 


+iG 9 »W 3 (T 013 - T m - T ml . 

Inspection of the right hand side of the equation above indicates that 
we need to calculate only 10 components of T afiv . Dispensing with terms of 
order (1 jr) 2 and higher orders, these components are given asymptotically 
by 


^201 

T;m 

Tl02 

T202 

TlQ3 

T;i0ii 

To 12 
T212 
71)13 
Til 3 


~ -d 0 1 + 


— —(d 0 l)sm8 + 


—cd 0 l — d 0 (cl + 2 Ml) — 2 dd 0 l + d 0 p 
sin 8 


ld 0 M + cd 0 l — d 0 (2Ml — cl + 2 dl) + d 0 p 


~ d Q l + 


- d 2 M + Mdol + e) 0 (2d + 2 dl + Ml) - ld 0 c - 2 ld 0 d 


—Idol H - rd 0 c cdoc Addod , 
(d 0 ) sin 0 
sin 8 


d 3 M + ld 0 c + Md 0 l + d 0 (M/ — 2 cl + 2 dl) — d 0 p 


~ — (/d 0 /+ rd 0 c — cd 0 c) sin 2 6*, 
doM 


d 2 l-M , 
d 3 M 


~ (d 3 Z) sin 6 — M sin 2 8 + l sin 8 cos 8 , 


( 22 ) 


No term of order 1/r in the expressions above will contribute to P^°h 
However, they will contribute to the momenta PT. The substitution of Eqs. 
m and (|22li into (|2Th yields 


s (°)oi 


M 


d 2 / + 


d 3 / cos 8 

~\~ Ij T~ 


sin 8 sin 8 


— i/d 0 / — ^ ld 0 l + cd 0 c + 2 dd 0 d 


(23) 


Assuming /(0) = /(0+27 t) and l sin 8 = 0 for 8 = 0 and 8 — n, it is easy to see 
that, under integration, the second term in the expression above vanishes, 
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r r d9d<f>wne(d 2 l + -^ + l^)=0. (24) 

Jo Jo V sm 9 sin 9 J 


Finally, substitution of Eq. 
the limit r —y oo, results in 


and e = r 2 sin 6 into definition (ITT]) . in 


r 2n rir 

p(°) = 4k dtj> dd sin 9 
Jo Jo 


M + d 0 F 


(25) 


where 


F = -i(l 2 + p) + + (26) 

The integral of M(u,9,(p) yields the mass aspect, but the total value of the 
gravitational energy depends also on the functions c(u, 9, 0) and d(u,9,(p), 
as one would a priori expect. The new term doF generalises the standard 
Bondi-Sachs energy. 

For finite values of the ordinary time t, the limit r — y oo corresponds to 
u —» —oo, and therefore in Eq. ()25]) we have P^°\u = —oo). However, this 
is not the general case. We may consider large but finite values of the radial 
coordinate r (r M and r ^ 8qF) such that Eq. (IT]) is verified. In this 
case, P < ‘°' ) depends on arbitrary and finite values of u. 

In equations (I22|) . we observe that T 2 02 and T 30 3 depend linearly on the 
radial coordinate r. This dependence could, in principle, lead to a divergent 
expression for P(°K However, the two contributions cancel each other, and 
the final expression turns out to be finite. 

In order to check the consistency of the result above, we calculated P ^ 
out of a slightly different set of tetrad fields that satisfies the same re¬ 
quirements that led to Eq. 031, namely, that the frame is adapted to a 
static observer and that the set of tetrad fields has the asymptotic form 
e at i — Va[i + (1/2 )h ai i in the limit r —>■ 00 . We considered e a/i given by 


e (o )n — ( — A —E, ~E, —G ), 

e (i> — (0, B\ sin 9 cos 0 + B 2 cos 9 coscp — B 3 sin# sin 0 , 
Cr cos 9 cos (p , 

—Dp sin 9 sin <p + D 2 cos 9 cos <p ), 
e ( 2 )/i = (0, B\ sin 9 sin 0 + B 2 cos 9 sin 0 + B 3 sin 9 cos 0, 
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Cr cos 9 sin <fi , 

Dir sin 9 cos </> + D 2 cos 9 sin </>), 

e( 3 )/i = (0, B\ cos 9 — B 2 sin 9, — C\r cos 9, — D- 2 sin 6*). (27) 

The values of A, B, ■ ■ ■ in the equation above are of course different from Eq. 
(USD, but eventually we obtain exactly as given by Eq. (125]) . The latter 
is, therefore, the gravitational energy that a static observer would measure 
at large distances from the source determined by M(u,9,(f)). Assuming the 
speed of light c = 1 as well as G = 1, the quantity 

Erad = ~r d(j) / d9 sin 9 (d 0 F) , (28) 

47T Jo Jo 

may be interpreted as the energy of the gravitational radiation. 

4.2 Gravitational momenta 

The evaluation of the gravitational momenta requires the calculation 
of £b) 01 — e w ^E^ 01 . The expressions of E^ 01 were already obtained in the 
evaluation of p(°>. These are the quantities that arise on the right hand 
side of Eq. (1211) . In the course of the calculations we find that the radial 
dependence of the tetrad components ew 2 and ep) 3 in Eq. (ITT)) requires the 
values of some components of T afiu of order 1/r, which are already presented 
in Eq. ()22l) . As we mentioned earlier, the contributions from the functions p 
and p cancel out in the calculations. After a number of simplifications and 
cancellations, we arrive at 


r 2n r n 

pW =4 k dcf) / d9 

Jo Jo 


sin 2 9 cos (f>(d 0 F) 


— ^sin 9 cos 9 cos <j>(2d 2 M + ld 0 M ) 


— sin + l sin 9doM) 


p ( 2 ) — 


p2 it 


sin 2 9 sm(p(d 0 F) 


4 k d(j) / d9 
Jo Jo 

+ - fsin 9 cos 9 sin (j)(2d 2 M + ld 0 M) 


+ cos 4>{2d 3 M + l sin 9d 0 M ) 
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(29) 


p( 3) 


4A) 


r 2rr 


deft / dd 


sind cos 6(d 0 F) 


■ sin 2 9(2d 2 M + ld 0 M ) 


The expressions above may be simplified by making integrations by parts in 
the angular variables. The expressions are further simplified by introducing 
the three-dimensional vectors, 


f * = (sin 6 cos 0, sin 9 sin 0, cos 9 ), 

9 l = (cos 9 cos 0, cos 9 sin 0, — sin 9 ), 

0* = (—sin0, cos0, 0). (30) 

We finally obtain 


P w = 4 k [ 2n d(j> r d9 sin 9 \(M + d 0 F)P 
Jo Jo L 

+-(ld„M)e‘ + 


(31) 


The expression above and given by Eq. (1251) constitute the gravita¬ 
tional energy-momentum of the Bondi-Sachs space-time in the realm of the 
TEGR. One important observation is the following. It is well known that the 
time derivative d$M may be written, in view of the field equations, in terms 
of the time derivatives of the functions c and d EM according to 

doM = —[(<9oc) 2 + (<%<i)"] + —dofdol + l cot 9 H— T — -'j . (32) 

2 V sm 60 

If we evaluate all integrals in P a in the limit r —>■ oo, we are actually taking 
the limit u — > — cxo. If we further assume that the news functions satisfy the 
initial conditions 


<9 0 c(—oo) = 0 , <9 0 <i(— 1 oo) = 0 , (33) 

then the total energy-momentum P a = (p(°),pW) given by (125|) and (|3Tj) 
reduces to the well known expression for the Bondi-Sachs energy-momentum. 
These initial conditions might be physically reasonable, but there is no justifi¬ 
cation for assuming them. We asserted above that we can evaluate (P®, pW) 
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over a surface S of integration sufficiently large from the source, such that 
Eq. (I23p is verified. This fact allows us to dispense with the initial conditions 
given by Eq. (j33jh 


5 News function in an axially symmetric space- 
time 

In this section we will discuss an application of Eq. (128]) by requiring the 
space-time geometry to be axially symmetric. The restriction to such sym¬ 
metry is achieved by enforcing d = 0 = l. The latter equations imply that 
the components of the space-time metric tensor do not depend on the vari¬ 
able 0, ensuring the axial symmetry of the configuration. This simplification 
allows the use a general form of the news function suggested by Papapetrou 
[20] . Halliday and Janis [21] and Hobill [22]. It is very difficult to obtain the 
exact form of the news function c{u, 9) because the latter must be consistent 
with the mass coefficient M{u , 6) such that Eq. (l32l) is verified, and the 
verification of the latter equation is non-trivial. 

Adopting the usual convention x = cos 9, the suggested form of the news 
function is written as [22] 



(34) 


where P n {x ) are the Legendre polynomials and a n {u) are general functions of 
the retarded time u, that reduce to constants when the source is not radiating. 
Of course, it is assumed that the summation in the expression above (as well 
as all summations in this section) converges. Hobill [22] specified a form of 
the time dependence of a(u ) and proceeded to obtain a particular expression 
for c(u,x). Below we will consider the more general form given by Eq. (134|) . 

The function F given by Eq. (1251) may be rewritten as 



(35) 


With the help of identities involving the Legendre polynomials, we write 
c(u, x) in the form 



n =2 


(36) 
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Substitution of Eq. (155|) into definition (125j) yields 


E r nrl . TV) 


f-1 r 


'-1 L 


■(1 — x 2 )(d x c) 2 + 4 xcd x c + 2c 2 (1 — 3x 2 ) 


dx. 


(37) 


The equation above may be further simplified by taking into account Eq. 
(I36jh We arrive at 


E r ad 0 d Q 

o 


f (1 - x 2 )(d x c) 2 dx - 6 a ndm f x 2 P 2 P^dx 

n,m 


(38) 


where P™ are the associated Legendre polynomials. Finally, using 


d x c = J2 


n=2 


p3 

n 


‘IxPl 


( 1-x 2 ) 1/2 (1-x 2 ) 


we obtain 


E, 


rad 


-kz 

* n 

4s«z 

° n,m 

4 a «z 


a 2 n (n + 2)!(16n 4 + 28n 3 + 20n 2 + lln - 105) 
(n - 2)!(2n + l)(2n - l)(2n + 3) 
q n q m (n + 2)![l + (-lf+ m ] 

(n — 2)! 

_ a n an+2(ri + 4)!_ 

(n - 2)! (2 n + 1)(2 n + 3)(2n + 5) ' 


(39) 


The Bondi mass m{u ) is defined as the usual integral over the mass aspect 


1 r 1 

m(u) = - j dx M(u, x). 

It has been demonstrated [20j !2U [22| that if one requires the total variation of 
the Bondi mass, due to gravitational radiation, to be equal to the total vari¬ 
ation of the mass aspect, then the summation in Eqs. (1381) and (1391) cannot 
have a finite number of terms, i.e., the summation cannot be truncated. 

The most interesting feature of Eq. (j39j) is the dependence of E rad on 
a quadratic combination of the coefficients a n (u). It might be possible to 
impose conditions on the coefficients a n (u) such that the summation simpli¬ 
fies considerably, but we have not attempted to carry out any simplification, 
which would require an analysis of the field equations. 
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6 Conclusions 


The Bondi-Sachs space-time is determined by three functions, M, c and 
d, which depend on The standard expression of the gravitational 

energy of the Bondi-Sachs space-time is the integral of the mass aspect 
M(u, 9, 0) in the angular variables, and restricting to the Bondi space-time, it 
is identified with m(u) given by Eq. (14Op . It is not clear why both expressions 
do not depend on the functions c and d. The expression that we obtained 
in the realm of the TEGR, Eq. ([25]) . does depend on M, c and d. We also 
obtained contributions of the functions c and d to the total gravitational 
momenta, as given by Eq. dBT]) . Altogether, Eqs. (j25|) and (IBT]) constitute 
the gravitational energy-momentum of the Bondi-Sachs space-time in the 
framework of the TEGR. 

The total energy-momentum of the Bondi-Sachs space-time, evaluated by 
means of the ADM definitions in (t, r, 9 , 0) coordinates, does depend on the 
functions c and d 112]. However, the ADM definitions are strictly valid only 
for asymptotically flat space-times, which is not the case of the Bondi-Sachs 
space-time. 

In view of the structure of Eq. (|25|) . we may identify the term that 
depends only on c and d, Eq. (|28|) . as the gravitational energy of radiation 
E radl since the functions c and d are non-vanishing at spacelike infinity and 
the news functions propagate over the whole three-dimensional space. The 
expression of E rad is clearly finite. 

By requiring d = 0 = l, we restrict the Bondi-Sachs space-time to the 
Bondi space-time, which is endowed with axial symmetry. In this space- 
time, a suggested form for the news function c(u,9 ) is given by Eq. (ITT) . 
The functions c and d yield the news functions, which are assumed as the 
radiating degrees of freedom of the gravitational field. From this point of 
view, it is interesting that the energy of gravitational radiation ([28]) results 
in a quadratic combination of the coefficients a n . 
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